Bulk viscosity and r-modes of neutron stars 



Debarati Chatterjee ( a ' b ) and Debades Bandyopadhyay( ab ) 

( a ) Theory Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Calcutta 
700 064, India 

( b ) Centre for Astroparticle Physics, Saha Institute of Nuclear Physics, 1/AF 
Bidhannagar, Calcutta 700 064, India 

Abstract. The bulk viscosity due to the non-leptonic process involving hyperons in 
K~ condensed matter is discussed here. We find that the bulk viscosity is modified 
in a superconducting phase. Further, we demonstrate how the exotic bulk viscosity 
coefficient influences r-modes of neutron stars which might be sources of detectable 
gravitational waves. 



1. Introduction 

Neutron stars have different classes of quasi normal modes (QNMs) [T] depending on 
restoring forces acting on a perturbed fluid element. Here we are interested in Coriolis 
restored fluid modes known as the inertial r-modes. The r-modes become unstable due 
to gravitational radiation. It might be responsible for regulating the spin of newly born 
neutron stars as well as old, accreting neutron stars. It was shown that the leptonic 
and non-leptonic bulk viscosities might effectively damp the instability in different 
temperature regimes. The bulk viscosity coefficient due to non-leptonic weak processes 
in neutron star matter involving hyperons [2j[3j|Il[5l[6l[7J[El[9] , antikaon condensate 
[TOl [IT] and quark matter [121 [131 [Hj was extensively calculated by many authors. 

In this article, we discuss the hyperon bulk viscosity due to the non-leptonic process 
n + p ^ p + A in K~ condensed matter and its role on damping the r-mode instability. 
In section 2, we describe the composition and equation of state (EoS) of neutron star 
matter. Results of hyperon bulk viscosity coefficients are explained in section 3. Section 
4 gives the summary. 



2. Composition and Equation of State 

Here we describe the composition and EoS of neutron star matter under going a first 
order phase transition from hadronic to K~ condensed matter. The constituents of 
charge neutral and /3-equilibrated matter in both phases are neutrons (n), protons (p), 
A hyperons, electrons, muons and also K~ mesons in the condensed phase. The baryon- 
baryon interaction is mediated by the exchange of cr, u and p, /o(975) (denoted hereafter 
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Figure 1. Bulk viscosity coefficient as a function of temperature with and without 
antikaon condensate for fixed values of baryon density (rib) and angular velocity (w). 



as a*) and 0(1020) mesons and described by the following Lagrangian density in Ref. 

Similarly we adopt the Lagrangian density for (anti)kaons in the minimal coupling 
scheme as given by (13, [TBI IH EDI ETj . The mixed phase is determined by the Gibbs 
phase rules and global charge neutrality and baryon number conservation |22j . 

Nucleon-meson coupling constants which are determined by reproducing nuclear 
matter saturation properties such as binding energy E/B = —16.3 MeV, baryon density 
no = 0.153 fm~ 3 , asymmetry energy coefficient a asy = 32.5 MeV, incompressibility 
K = 300 MeV and effective nucleon mass m* N /m^ = 0.70, are taken from Ref.[23j. 
Further kaon-meson and hyperon-meson coupling constants are obtained from Ref. 
[TTl fT5] . In this calculation, we adopt the value of antikaon optical potential depth 
at normal nuclear matter density as Ux(no) = — 160 MeV. 

When there is no first order phase transition from hadronic to K~ condensed 
matter, A hyperons appear at density 2.25no- On the other hand, the early onset 
of K~ condensate in a first order antikaon condensation, delays the appearance of A 
hyperons. In our calculation, K~ condensation begins at 2.23n and the mixed phase 
ends at 4.1n whereas A hyperons appear at 2.51n . Further we calculate energy density 
and pressure corresponding to above mentioned compositions. 
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3. Hyperon bulk viscosity coefficient and suppression of r-modes 



The real part of bulk viscosity coefficient due to the non-leptonic process n + p ^ p + A 
is given by [U [21] , 

1 + {ujrf 

The difference of infinite and zero frequency adiabatic indices is, 

_ _ _r^_dP_dXn , * 

^° P <9n„ drib ' 

where x n = ^ is the neutron fraction in the equilibrium state. The relaxation time (r) 
for the non-leptonic process in the j-th (= hadronic / K~ condensed) phase is given by 

2 

where p\ is the Fermi momentum for A hyperons and < \M\\ > is the angle averaged 
matrix element squared in the corresponding phase. 

Temperature dependence of hyperon bulk viscosities due to the non-leptonic process 
n + p f± p + A with and without K~ condensate are shown at rib = 0.45/m -3 and 
uj = 10 4 s _1 in Figure 1. The solid line represents the hyperon bulk viscosity in the 
absence of antikaon condensate. The dashed and dashed-dotted lines correspond to the 
hyperon bulk viscosities in the hadronic and antikaon condensed phases of the mixed 
phase in a first order hadronic to K~ condensed phase transition [If J. The hyperon bulk 
viscosity in the absence of K~ condensate is suppressed compared with hyperon bulk 
viscosities in the presence of antikaon condensate. Further we note that the hyperon bulk 
viscosity coefficient in the condensed phase is smaller than that of the hadronic phase 
above 3 xf0 8 K. This is also the temperature where the inversion of the temperature 
dependence of hyperon bulk viscosity in K~ condensate happens. This feature is also 
found in the case without the condensate. 

We investigate the r-modes of a neutron star having gravitational mass f.6 M 
corresponding to central baryon density 3.5no and rotating at an angular velocity 2652 

We calculate damping time scales corresponding to hyperon bulk viscosity (r BV ), 
modified Urea bulk viscosity involving nucleons (tmu), shear viscosity (tsv) and the 
growth time scale of gravitational radiation (tq) using the energy density and bulk 
viscosity profiles [TT]. We compute the critical angular velocities (flc) as a function 
temperature solving — - — h — + — H — — = . Various damping time scales are 
exhibited as a function of critical angular velocity in Figure 2. It is noted that the 
r-mode instability of the neutron star is damped when the growth time scale due to 
gravitational radiation is comparable with damping time scales tbv and Tmu- We find 
that the hyperon bulk viscosity damps the instability at T = 4 x f 9 and below whereas 
it is the modified Urea bulk viscosity which effectively suppresses the r-mode instability 
above this temperature. 
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Figure 2. Different damping time scales and the growth time scale due to gravitational 
radiation are shown as a function critical angular velocity (f2c)- 



4. Summary 

We have studied hyperon bulk viscosity in the presence of K~ condensate. We find the 
inversion of the temperature dependence of hyperon bulk viscosity in our calculation. 
Though the hyperon bulk viscosity is suppressed in K~ condensed matter than that in 
the hadronic phase, it still effectively damps the r-mode instability. 
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